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Charged lepton flavor violating processes are unobservable in the standard model, but they are 
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In the original formulation of the standard model (SM) 
in which neutrinos are massless, lepton flavor is an ac- 
cidentally conserved quantum number. The extension 
of the SM to include neutrino masses introduces os- 
cillations between the neutrino flavors, which violate 
this conservation law. However, SM processes involv- 



ing charged lepton flavor violation (CLFV) remain un- 
observable because they are suppressed by the quantity 
{Aml/M^f < 10-48 Here Ami is the difference 

between the squared masses of neutrinos of different fla- 
vor and Mw is the charged weak vector boson mass. 
Hence CLFV represents an unambiguous signature of 
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new physics (NP) 0-0|- Many extensions to the SM, in- 
cluding supersymmetry and models with leptoquarks or 
compositeness, predict an enhancement in the rates for 
these processes at levels close to experimental sensitiv- 
ity. There have been considerable efforts in searches for 
CLFV in decays of particles such as /i and r leptons and 
B and K mesons, but CLFV in the T sector remains rel- 
atively unexplored Q. By using unitarity considerations, 
limits on CLFV r branching fractions [9| have been used 
to place indirect limits on CLFV T branching fractions at 
the 0(10"^) level [Ifl]. In this letter we describe a search 
for CLFV T decays, which is a thousand times more sen- 
sitive than these indirect limits, using data collected with 
the BaBAR detector at the PEP-II B factory at SLAC Na- 
tional Accelerator Laboratory. Since these decays are in 
general mediated by new particles produced off-shell in 
loops, their measurement probes mass scales up to the 
TeV scale, far exceeding the e+e~ center-of-mass (CM) 
cohision energy = Mr(nS) = O{10 GeV) There- 
fore, this analysis provides a NP probe which is comple- 
mentary to direct searches ongoing at the Tevatron and 
to be performed at the Large Hadron Collider. 

Assuming that the partial widths for CLFV T decays 
are comparable at the T{2S), T{3S) and T(4S') reso- 
nances, the branching fractions for rare decays of the 
narrow T(jiS) resonances (henceforth n = 2,3) are en- 
hanced by approximately Tr(4s)/T'T(nS) — 0(10"^) 
with respect to those of the T{AS). We search for the 
CLFV decays r{nS) e^r^ and r{nS) ->■ h^t^ [H, 
while the decay T{nS) — >■ e'^fi^ is constraine d by uni- 
tarity considerations to be less than 0(10"*) [10|. No 
signal is expected in data collected at the T{4S) since 
the CLFV branching fractions are strongly suppressed, 
or in data collected away from the T resonances, since 
this data contains very few T decays. We search for 
CLFV in a sample of (98.6 ± 0.9) x 10*^ T{2S) decays 
and (116.7 ± 1.2) x 10^ T{3S) decays corresponding to 
integrated luminosities of 13.6 fb""'^ and 26.8 fb~^, respec- 
tively. Data collected at the Y{AS) after the upgrade of 
the muon detector system (77.7 fb~^) and data collected 
30 MeV below the T{2S) and T{3S) resonances (off-peak 
data corresponding to 2.6 fb"^ and 1.3 fb^^, respectively) 
constitute control samples that are used to validate the 
fit procedure. An additional data control sample col- 
lected at the T(3S') resonance (1.2 fb^^) is used in a 
preliminary unblinded analysis to validate the analysis 
procedure and to ensure agreement between data and 
events simulated using Monte Carlo (MC) techniques. 
Simulated background processes consisting of continuum 
QED events [3, [11] and generic T{nS) decays, as well 
as signal T{nS) — >■ ("^t^ {£ = e, fi) decays [16[, are pro- 
duced and analyzed to optimize the fit procedure. The 
Geant4 |17i] software is used to simulate the interac- 
tions of particles traversing the BaBAR detector, which is 
described in detail elsewhere fi^. 

The signature for T{nS) £^t^ events consists of ex- 



actly two oppositely charged particles: a primary lepton, 
an electron (muon) for the T{nS) — )• e^r^ {T{nS) — )■ 
/i^r^) search, with momentum close to the beam energy 
Eb — -\/s/2, and a secondary charged lepton or charged 
pion from the r decay. Here and in the following all 
quantities are defined in the CM frame unless otherwise 
specified. If the r decays to leptons, we require that the 
primary and r-daughter leptons are of different flavor. If 
the r decays to hadrons, we require one or two additional 
neutral pions from this decay. These requirements on the 
identified particle types are necessary to suppress Bhabha 
and /x-pair backgrounds. Thus we define four signal chan- 
nels, consisting of leptonic and hadronic r decay modes 
for the T{nS) — ?> e'^r^ and T{nS) — > fi^r^ searches, 
hereafter referred to as the leptonic and hadronic er and 
fiT channels. The main source of background to our 
events comes from r-pair production, for which the final 
state particles are the same as for the signal. There is a 
background contribution to the er channels from Bhabha 
events in which one of the electrons is misidcntified, and 
to the fiT channels from /i-pair events in which one of 
the muons is misidcntified or decays in flight, or an elec- 
tron is generated in a material interaction. An additional 
background consisting of events with multiple pions and 
possible additional photons ('7r-hadron background'), in 
which a charged pion is misidentified as a lepton and the 
remaining particles pass the selection criteria for the r 
decay products, contributes to the hadronic er and /ir 
channels. 

In order to reduce background, we first apply require- 
ments common to all the decay modes and then a channel 
specific selection. All events are required to have ex- 
actly two tracks of opposite charge, both consistent with 
originating from the primary interaction point and with 
opening angle greater than 90° . To suppress Bhabha and 
/i-pair backgrounds, we require that Myis/^/s < 0.95, 
where M^s is the invariant mass of the sum of the 4- 
vectors of the two charged particles and of all photon 
candidates in the event. To ensure that the missing mo- 
mentum is not pointing toward the holes in the detec- 
tor near the beamhne, we require that cos(6'^^gg) < 0.9 
and cos(Cfi.) > -0.9, where Ol^l, (CffJ is the po- 
lar angle of the missing momentum in the lab (CM) 
frame. To suppress two-photon processes, we require that 
(Pi + P2)±/{Vs - \Pi\ ~ > 0.2, where pi and P2 
are the momenta of the two charged particles and _L indi- 
cates the transverse component with respect to the beam 
axis. 

Particle identification is performed using a multivari- 
ate analysis [l^ which uses measurements from all of the 
detector subsystems. An electron selector and a muon 
veto, combined with the requirement that the particle 
falls within the angular acceptance of the electromag- 
netic calorimeter (EMC), are used to identify electrons. 
A muon selector and an electron veto are used to iden- 
tify muons, while a charged pion selector, an electron 
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veto and a muon veto are used to identify charged pions. 
The particle misidentification efficiencies are 0(10"^) 
ifi e), O{10'^) (e ^ ^) and O{10-^) {£ tt), where 
X — >■ 2/ indicates that a particle of type x is misidenti- 
fied as a particle of type y. A photon candidate must 
deposit at least 50 MeV in the EMC and have a shower 
profile consistent with that expected from an electromag- 
netic shower. All pairs of photons with an invariant mass 
between 0.11 GeV and 0.16 GeV are selected as neutral 
pion candidates. 

The channel specific selection classifies events into one 
of the four signal channels. The momentum of the pri- 
mary lepton normalized to the beam energy is required 
to satisfy x = \pi\/Eb > 0.75. For the hadronic r-decay 
channels, the momentum of the r-daughter charged pion 
is required to satisfy \p2\/Eb < 0.8. Since these r de- 
cays to hadronic final states are dominated by the de- 
cays — > p^i^T and — !■ afi/r, the masses of the 
7r^7r° and tt^tt^tt^ systems are required to be consistent 



leptonicT(3SH ex (x /ndf=52.4/49) 



with the masses 



0.77 GeV and 



= 1.26 GeV, 



respectively, where the requirement on the n^n^n'^ sys- 
tem mass is included only if there are two neutral pi- 
ons in the event. In order to suppress Bhabha events 
in which an electron is misidentified as a muon, for the 
leptonic er channel the r-daughter muon is required to 
penetrate deeply into the muon detector. In order to 
suppress /i-pair events in which the tracks are back-to- 
back, for the leptonic /xr channel the tracks are required 
to satisfy A(j) < 172°, where A(j> is the difference between 
the track azimuthal angles. After including all selection 
requirements, typical signal efficiencies determined from 
MC are (4 — 6)% [23|, including the r decay branching 
fractions. The typical number of events passing the se- 
lection criteria is (10 — 15) x lO'^ for T{2S) data and 
(20 - 30) X 10^ for T{3S) data, depending on the signal 
channel. These yields are consistent with background 
expectations from MC simulations. 

After selection an unbinned extended maximum likeli- 
hood fit is performed to the distribution of the discrim- 
inant variable x. The signal peaks at a; w 0.97, while 
the T-pair background x distribution is smooth and ap- 
proaches zero as X ^ xmax, where xmax ~ 0.97 is the 
effective kinematic endpoint for the lepton momentum 



in the decay 



-i^iVt, boosted into the T{nS) rest- 



frame. The X distributions for the Bhabha///-pair back- 
grounds have a peaking component near x — 1, about 
(2.5-3)(Ta; above the signal peak, where ax ~ 0.01 denotes 
the detector x resolution. The x distribution for the tt- 
hadron background is smooth and falls off sharply near 
X = Xmax- Probability density functions (PDFs) for 
signal, T-pair, Bhabha//x-pair and 7r-hadron backgrounds 
are determined as discussed below, and a PDF consist- 
ing of the sum of these components weighted by their 
yields is fitted to the data for each signal channel, with 
the yields of the components allowed to vary in the fit. 
The PDFs for signal and Bhabha//x-pair backgrounds 
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FIG. 1: Maximum likelihood fit results for the leptonic er 
channel in T{3S) data. The red dotted line represents the 
signal PDF, the green dashed line represents the sum of all 
background PDFs and the solid blue line represents the sum 
of these components. The inset shows a close-up of the region 
0.95 < a; < 1.02. The top plot shows the normalized residuals 
(data - fit) /(Jdata (puU). 

are extracted from fits to the x distributions of MC 
events. The signal is modeled by a modified Gaussian 
with low- and high-energy tails, hereafter referred to as 
a double Crystal Ball 2l| function, which peaks near 
X = 0.97. The Bhabha and //-pair backgrounds have a 
threshold component truncating near x = 1, which is 
modeled by an ARGUS distribution [g^l, and a peaking 
component near x = 1, which is modeled by a Gaus- 
sian function. The 7r-hadron PDF is determined from 
data by modifying the selection to require that the pri- 
mary lepton is instead identified as a charged pion. The 
resulting binned x distribution is scaled by the proba- 
bility for pions to be misreconstructed as charged lep- 
tons, as measured in data, to yield a binned PDF for the 
TT-hadron background. The yield of this component is 
fixed in the maximum likelihood fit and an uncertainty 
of 10% is assessed. The r-pair background is modeled by 
the convolution of a polynomial, which vanishes above 
the kinematic endpoint Xmax , and a detector resolution 
function. The detector resolution function is modeled by 
a double Crystal Ball function whose shape is extracted 
from T-pair MC events. Since the signal peaks in the 
region near the kinematic endpoint of the r-pair back- 
ground X distribution, the signal yield depends strongly 
on X MAX , which must therefore be extracted from data. 
The value of this parameter is extracted from fits to the 
T{4S) data control sample and corrected for differences 
in the decay kinematics at the T(4S') vs. T{nS) res- 
onances. The polynomial shape parameters, which are 
not strongly correlated with the signal yield, are allowed 
to vary in the fits to T{nS) data. 

To validate the fit procedure, we perform fits to data 
control samples in order to verify that signal yields con- 
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sistent with zero are obtained. The T(4S') data is divided 
into samples that are chosen to be comparable in size to 
the T{2S) and T{3S) data samples. The off-peak data 
and the 1.2 fb~^ of T(3S') data constitute additional data 
control samples. Results consistent with zero signal yield 
are obtained for all signal channels in these data control 
samples. 

The branching fraction B is calculated from the 
extracted signal yield Nsig according to B = 
NsigI {^SIg X Nr{nS))i where esiG is the signal selection 
efficiency and A^r(nS) is the number of collected T{nS) 
decays. The dominant systematic uncertainties in the 
signal yields, which arise from uncertainties in the PDF 
shapes, are determined by varying the shape parameters 
while taking into account the correlations between them. 
This uncertainty is 3-10 events depending on the signal 
channel, and the largest contribution is due to the un- 
certainty in the kinematic endpoint parameter Xmax- 
To assess the uncertainty in the signal efficiency, we take 
the relative difference between the yields for data and MC 
events from a portion of the sideband of the x distribution 
defined by 0.8 < x < 0.9, which is dominated by r-pair 
events. This difference is due to particle identification, 
tracking, trigger, and kinematic selection efficiency un- 
certainties. There is an additional statistical uncertainty 
in the signal efficiency, as well as an uncertainty arising 
from the uncertainties in the t branching fractions [l3 |. 
The total signal efficiency uncertainties are (2 — 4)%, de- 
pending on the signal channel. The uncertainty on the 
number of collected T{nS) decays is approximately 1%. 

To assess the possible bias in the fit procedure, sev- 
eral hundred simulated experiments are produced with 
the generated signal yield fixed to the larger of the value 
extracted by the fit to T{nS) data, or zero. The bias 
is consistent with zero within the uncertainty of 0.2-0.7 
events, depending on the signal channel. There is also 
an uncertainty resulting from a correction in the signal 
yield which is performed to compensate for primary lep- 
tons whose momentum is poorly measured. These par- 
ticles populate a broad momentum range and some fall 
in the signal region defined as the interval within ±1.5f7 
of the signal peak. The number of these events is es- 
timated using T-pair MC simulation and corrected us- 
ing data and MC Bhabha and ^-pair control samples. 
The expected contributions are subtracted from the sig- 
nal yields extracted by the fit and an uncertainty of 100% 
times the correction is assessed. The corrections are ap- 
proximately 3 events (5 events) for the T{2S) — > /i^r^ 
{T{3S) — ^ /i^T^) channels and less than 1 event for the 
T{nS) -> e^T^ channels. 

The maximum likelihood fit results for a sample chan- 
nel are displayed in Fig.^l and the fit results for all 
channels are available in [20j. After including statistical 
and systematic uncertainties, the extracted signal yields 
for all channels are consistent with zero within ±1.8(T. 
We conclude that no statistically significant signal is ob- 



TABLE I: Branching fractions and 90% CL ULs for signal 
decays. The first error is statistical and the second is system- 
atic. 



B (IQ-^) UL (10-^) 



B{T{2S) - 


e±r=F) 


n « + l-5+0.5 
"•"-1.4-0.6 


< 3.2 


B{r{2S) - 




r, 9 + 1.5 + 1.0 
"•^-1.3-1.2 


< 3.3 


B(r{3S) - 


e±rT) 


-'-•"-1.4-0.7 


< 4.2 


B(r(3S') - 




n + 1.5+1.4 


< 3.1 



served and determine 90% confidence level (CL) upper 
limits (UL) using a Bayesian technique, in which the 
prior likelihood is uniform in B and assumes that B > 0. 
The resulting ULs, summarized in Table HI are 0(10"^) 
and represent the first constraints on B{T[nS) — ?> e^r^). 
These results improve the sensitivity by factors of 3.7 
and 5.5, respectively, with respect to the previous ULs 
on B{T{2S) n'^T^) and B{r{3S) -> fi'^T^) 

Our results can be used to constrain NP using effective 
field theory. The CLFV T{nS) decays may be parame- 
terized as an effective bM^r^ 4-fermion interaction given 
by [81] 

AC = ^iiT,r)ih^b), (1) 

where F^ is a vector or an axial current or their combina- 
tion, air and Air are the NP coupling constant and mass 
scale, respectively. This allows the following relation to 
be derived [Uliil: 

ajr _ B{Y{nS) ^ (.^t^) 2gga^ 
Aj^ ~ S(T(n5) -> £+£-) (Mr(„s))*' 

Here = —1/3 is the charge of the h quark, a = 
a{MY{nS)) is tfie fine structure constant evaluated at the 
T{nS) mass, and we take the dilepton branching frac- 
tion B{T{nS) -> I'^a^) from the average of the T(nS) 
dielectron and dimuon branching fractions [T2|. Us- 
ing these values and taking into account the 7% uncer- 
tainty in B{T{nS) — > £^£^), we determine the likelihood 
as a function of the quantity a'j^/Aj^ and extract the 
90% CL UL using the same Bayesian method as above. 
We use these results to exclude regions of the Air vs. 
air parameter spaces as shown in Fig. [2] Assuming 
aer = ctfxr = 1, these results translate to the 90% CL 
lower limits Aer > 1.6 TeV and yl^r > 1.7 TeV on the 
mass scale of NP contributing to CLFV T{nS) decays, 
which improve upon the previous lower limit on A^^r 18[. 
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FIG. 2: Excluded regions of effective field theory parameter 
spaces of mass scale Atr versus coupling constant ae-r. The 
dotted blue line is derived from T{2S) results only, the dashed 
red line is derived from T{3S) results only, and the solid black 
line indicates the combined results. The yellow shaded regions 
are excluded at 90% CL. 
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APPENDIX: EPAPS MATERIAL 



The following includes supplementary material for the Electronic Physics Auxiliary Publication Service. 



TABLE IL Signal efficiencies esic and signal yields Nsig extracted by the maximum likelihood fit for the four signal channels 
in T(25) and T{3S) data. The first error is statistical and the second is systematic. 



leptonic er hadronic er leptonic fir hadronic /ir 

r(2S) esiG (%) 5.10 ±0.03 ±0.10 5.53 ± 0.03 ± 0.16 4.22 ± 0.03 ± 0.11 6.08 ± 0.03 ± 0.12 

Nsig -8 ± 8 ± 3 19 ± 12 ± 4 -11±8±5 17 ±13 ±7 

T{3S) esiG (%) 5.24 ±0.05 ±0.17 5.56 ± 0.05 ± 0.12 4.20 ± 0.05 ± 0.06 6.07 ± 0.06 ± 0.10 

Nsig 24 ±13 ±4 -6 ± 13 ± 6 -11 ±10 ±9 4 ±16 ±11 



leptonicr(2S)-> ei (x^/ndf=53.2/49) 




0.75 0.8 0.85 0.9 0.95 1 



hadronicT(2S)^ ei (x^/ndf=31.6/49) 
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leptonicT(3S)^ ex (x^/ndf=52.4/49) 
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hadronicr(3S)-» ex (x^/ndf=40.6/49) 
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leplonlcr(2S)^ HI (x^/ndf=48.3/49) 



1 1^^^^^^^ 




-2 
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leptonic r(3S)^ HX (x /ndf=35.7/49) 




hadronicT(2S)^ (ix (x /ndf=32.9/49) 



hadronic r(3S)^ |ix (x /ndf=48.9/49) 
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FIG. 3: Maximum likelihood fit results for the four signal channels in T{2S) data (left) and T(3S) data (right). The red dotted 
line indicates the signal PDF, the green dashed line indicates the sum of all background PDFs and the solid blue line indicates 
the sum of these components. The inset shows a close-up of the region 0.95 < 2: < 1.02. The pull denotes the difference in 
each bin between the data and total PDF, normalized to the statistical uncertainty in data. The corrected signal yield Nsig 
is displayed with statistical and systematic errors. 



